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Abstract 

We calculate cross-sections for neutral current deep inelastic scattering at HERA with 
photon tagging. Both, the exact lowest-order cross-section and a leading logarithmic 
approximation of next-order corrections are calculated. The latter amounts to less than 
20% in a large kinematical range. The integrations are performed numerically but without 
relying on Monte Carlo methods. 
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1 Introduction 



We calculate cross-section predictions for the reaction 

eih) + pipi) e{k2) + X(p2) + l{k), (1.1) 

where the photon is assumed to be observed (tagged). The calculated cross-section will 
be differential in energy E'^ and angle 9'^ of the final state electron and in the energy EZ^"'^ 
of the photon. Equivalently, the usual scaling variables x, y together with E^^^ may be 
chosen^. The angular range covered by the photon tagging device is taken into account by 
kinematical cuts in the course of integrations. In the numerical results we will follow the 
experimental set-up of the HI collaboration with the photon tagging device located in the 
forward direction with respect to the electron beamQ. The lowest-order cross-section is 
determined without approximations. Since one expects large photonic corrections in some 
kinematical ranges, we tried to estimate them with a leading logarithmic approximation. 
What we do is certainly not a rigorous approach and cannot replace a complete second- 
order calculation. In this connection the neglection of effects due to photon interferences 
has to be mentioned. They are of next-to-leading order, but may be non-negligible for 
a forward tagging device. For the time being, we hope to cover the most substantial 
corrections. 

In Section |, the lowest-order cross-section is calculated. Section ^ contains a discus- 
sion of the inclusive cross-section. The leading-logarithmic corrections are calculated in 
Section || and numerical results are presented in Section |. Some technical details may 
be found in the Appendix. 

2 The Lowest-Order Cross-Section 

In this section we derive a set of formulae, which allow the lowest-order description of 
deep-inelastic cross-sections with tagged photons, where the latter may have arbitrary 
kinematics. Obviously, the treatment of arbitrary cuts in energy and production angles 
is possible only with numerical methods. Thus, we decided to tackle the problem not 
with the semi-analytical approach as developed in our previous papers. Instead, we will 
numerically integrate the squared matrix element over the phase space imposing the 
appropriate cuts explicitly. For this aim, we use natural variables in the HERA system. 
Besides the four-momenta of initial state electron and proton, this are variables, which 
are directly measured by the apparatus: angles and energies of final state electron, E^,9'^, 
and of the photon, E^,9-y,ip-y. The following representations for the 4- vectors pi,ki,k2,k 
have been chosen (see also Figure |l]): 

pi = (Ep,0,0,-pp), (2.1) 

fcl = (Ee,0,0,Pe), (2.2) 



^ We note that the usual definition of x and y in terms of leptonic variables is used. They are not 
recalculated using the reduced electron beam energy after emission of the observed photon. 

This direction is called backward direction [with respect to the proton beam] in the HI collaboration. 
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k2 = {E'^,p'^sme'^,0,p'^cose'^), (2.3) 

k — E.y{l,sm9^ cosip^,sm9^smip^,cos9^), (2.4) 

P2 = pi + ki - k2 - k. (2.5) 
In these variables, the invariants s, Qf, yi are: 

s = (ki+pif = {E, + Epy-{p,-ppy = S + M'' + m\ (2.6) 

S = 2pi/Ci = 2(Ee^p+PePp), (2.7) 

yi = ^^i^^^^ = l-2(EX+PpP>sO/'5, (2.8) 

= -{ki-k2f = 2{EeE'^-p^p'^cose'j-2m\ (2.9) 

With m and M we denote the electron and proton masses, respectively. We will also use 
the TT production threshold, 

M = M + m^, (2.10) 
which is the lowest mass of the inelastic final hadronic system. 
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The angular variables are unlimited, 

-1< cos^l < +1, 



:2.111 



-1 < cos 6*^ < +1, 



(2.12) 



< < 27r, (2.13) 

while the allowed ranges of electron and photon energies have to be determined: 



(2.14) 



0< E, <E^^''{E'^,e'^,e„^,). 



The upper limits are derived in Appendix ^ 



max 



with 



So 



{Ee + Ep){EeEp + pePp + m^- AMV2) + (pe - Pp) COS e'jLe 



Le = iE,Ep + p,pp + m' - -AM'Y -m'Se 



Se = {E, + Epf -{p,-ppf cos" e',, 



(2.15) 
, (2.16) 

(2.17) 
(2.18) 



and 



^7 



S + 2m' + KP - A?2 -2{E, + Ep)E[ + 2{p,- pp)p^ cos 6'^ 



2 [Ef. + Ep- E'^ + (pg cos 9'^ - [pe - Pp\) cos 6^ + p'^ sin 9'^ sin 6^ cos (f^] 
The phase space is defined as follows: 



dr 



dh dk dp2 J , ,2 c4 



2k^ 2A;0 2p! 



±dM^S\ki+pi-k2~k~p2), 



with 



Ml 



(2.19) 
(2.20) 

(2.21) 
(2.22) 



(pi + ki-k2~ kf. 

We show in appendix ^ that the doubly-differential phase space in terms of natural 
variables takes the following form: 



dV 



TT 



^max 



f-27r ^ + 1 

, duD^ / dcos6^ I E^dE^. 
p'^EJcosO^ 2 Jo ^ 'J^i ^ ^ 

The totally differential cross-section in natural variables is: 

d C^brem 



4q;3 , X ^ 



(2.23) 



(2.24) 



dcos6^dE^dcos6^d(p^dEy nS Qf^ -^^ 

The functions JFj are generalized structure functions. Explicit expressions for them may 
be found in Equations (1.2.18) 0. 

Here and henceforth, we mark equations taken from reference 111] with I. 
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The totally differential 'radiator functions' Si for the radiative process are, with ac- 
count of minor improvements and modifications, taken from Equations (I.2.40)-(I.2.42): 



Si 
S2 



1 n Qt + Qi 



(2.25) 



S\Ayi, + y,)-M'Ql , S\l-y,)-M'Qt 



z- 



^2 



S. 



1 

Vh 



Z1Z2 

22/;i +yh , '2-yh 



^2 
^2 



-S^ 







-2M2| 




^2 / 





(2.26) 



2^1 2^2 



-Ql ' ^" 



1 



Zl Z2 



Vh 



Zl Z2 



(2.27) 



Above, we retained m where needed and M throughout. The following abbreviations and 
definitions have been used: 



Vh 
Vh 

Zl 



1 



IE 

yi ^(^^p + PpCos6'-^), 



2kik = 2E^ {Ee — PeCos 6^) , 



(2.28) 
(2.29) 

(2.30) 
(2.31) 



Z2 



2k2k = zi + Qf-Ql 



2E^ 



E'^ — Pg (^cos 6'gCos 6^ + sin 6'gSin 6'^cos 



^1 



(2.32) 



Ql = -{P2~Pif 

= Qf — 2Ej {pe — PgCos 6'g)cos 9^ — PgSin 6'gSin 6'^cos (fj + E'^ — Ee . (2.33) 

The definitions reflect the radiative kinematics: the difference Ql — Qf as well as zi and 
Z2 are proportional to E^. In Appendix ^ the same is shown for — M^. 

From the above expressions, any radiative cross-section, if measured in natural vari- 
ables, may be calculated by simply imposing the corresponding cuts when integrat- 

ing ( ra . 

The lowest-order radiative cross-section becomes: 



dcos O'^dE'^dE^ nS 



Pg [ dcos6y I dip^^^SiTi 

i=i 



(2.34) 
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with reasonably defined boundaries for tlie angular integrations. In case of a forward 
tagging device as being used by the HI and ZEUS collaborations: 



cos 6*^ G 



-Ml, 1 



ip^ e (0,27r). 



(2.35) 



Here, A9^ describes the opening angle of the tagging device^. 
For applications, it may be useful finally to change variables: 



dxidyidE^ 



J 



d^CThr 



dcos OAEAE^ 



(2.36) 



with Ql = xiUiS. The Jacobean is easily determined from ( p^) and ( p^ 



J 



d{cos 6'^, E[ 



d{xi,yi] 



yiS 

2Pe 



(2.37) 



We show the lowest-order cross-sections ( p.36| ) in Figures ^ and |^ for E^ = 5, 10 GeV 



and A9.y = 0.5 mrad. The updated Fortran program HECTOR [0 with flag settings ISSE=1 
and ISCH=0 (selecting structure functions CTEQ3(L0)) is used for all numerical results. 
As is seen in the figures, the minimal value of y depends on both x and i?™ = E.y. 
In addition, it depends on the angular cuts. For a forward tagging device, its absolute 
minimum may be approximated by 



E^ 
E. ■ 



(2.38) 



The cross-section resembles the typical behavior of inclusive deep inelastic scattering. 



3 The Inclusive Cross-Section 

In the preceding Section, we made the assumption that the energy of the tagged photon 
will be observed. Events with vanishing photon energy were automatically excluded and 
thus all the problems related to the infrared singularity. In this Section, we will treat the 
photon energy as inclusive variable, i.e. determine the doubly-differential cross-section 
d'^a/dcos6^dE^ as the sum of Born cross-section, vertex correction, and real bremsstrah- 
lung contributions: 

<5.ert) + / dE^ ^^ ,_a',T.' - (3-1) 



d COS 9'^dE'^ d COS 9'^dE'^ ^^'^ J ^ dE^d cos 6'^dE'^ 

The purpose is two-fold. We want to have the opportunity to perform a numerical compar- 
ison of the present, deterministic calculations with those performed in the semi-analytical 
approach; and the latter are doubly-differential. Further, in the next Section we estimate 



^ Strictly speaking, one should take into account the actual form of the tagger (being e.g. rectangular 
rather than circular); we leave this as an exercise to the reader. 



6 




Figure 2: Lowest-order cross-section for HERA kinematics with — 5 GeV and A^-y = 
0.5 mrad 
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higher-order QED corrections using the leading-logarithmic approximation for them. In 
doing so, we will have to make use of ( |2.34|) without a restriction on the lower limit of E^. 



We will consider ( |2.34|) without cuts on the photonic angles 



dcos O^dE^dE^ nS 



rl r2-K 

p'^ / dcos9j / J VS'iJ'i, (3.2) 

J-i Jo Qt ^ 



and derive the lowest-order QED corrections d^a^^^^/ dxidyi to the cross-section 

dxidyi dxidyi dxidyi 

Then, d'^a^j^^^ has to be compared to the semi-analytical result, see e.g. Equations (1.4.40) 
and (1.5.3). For this aim, we have to multiply ( p.2| ) by the Jacobean ( [A.8| ), then to 
integrate over E^, and therefore to regularize the infrared divergence. From (|3.1| ), the 
^^"^QED t)e calculated now: 



d COS 6'^dE'^ 7-1 Jo Jo yd cos 6'^dE'^d cos O^dip^dE^ 

1 d^aB. 



jorn 

An^dcosO'JE'^ 



1 + 6,,{E,) 6{E,)}, (3.4) 



with a cutoff E^ on the minimum photon energy for the numerical integration. The 



correction S^-^ will be discussed below. Further, it is 



An^ = / dcos6j / d(p^. (3.5) 



In ( |3.4|) , the angular integrations are unrestricted, AQ^ = An and (i^o"brem is a completely 
differential, infrared-finite bremsstrahlung cross-section and may be treated numerically 
as it stands. It is the analogue to Equation (1.5.3) without the subtracted terms. 

In fact, we have to follow a different strategy for the solution of the soft-photon 
problem. In the semi-analytical approach, we performed a subtraction of the infrared 
divergent part of the integrand thus making the numerical integration infrared finite in 
the limit Ej ^ 0. This subtraction was compensated by an analytical integration of 
the subtracted part. We specially had to choose a reference frame where the limits of 
photon energy variation are independent of the angular variables, this way considerably 
simplifying the integrations. This was possible since the net cross-section was defined in a 
covariant way. Here, the situation is quite different. We have to use the above defined set 
of natural variables in order to finally be allowed to apply the cuts. As was discussed, the 
photon energy varies in limits which depend in a complicated way on the photon angles. 
It is important for the calculation that one may choose E^ so small that it is well inside a 
region where E^ is unrestricted by the photon angles. That this region exists for arbitrary 
values of 6^, E^ may be seen in Figure § where the allowed region for photon energies is 
shown below the shaded surface as function of these angles. At cos 6'^ = 0, the minimal 
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£;max -g reached. If is smaller, there is no influence of the photon angles on the photon 
energy boundary. We have to choose: 

g + 2m^ + - - 2{E, + Ep)E'^ + 2{p, - pp)p, cos O'^ 

2[i?, + E,-E; + (p;cos^;-[p,-p,])] • 

We split the correction into pieces as it was done in (1.4.40): 

" ^QED « r " (JBorn , O C^soft , " "^hard , " ^Tr 



d COS 9'^dE'^ TT d cos 9'^dE'^ d cos 9'^dE'^ d cos 9'^dE'^ d cos 9'^dE'^ 



- r-SvalsJrj^r^^ + rjT^^, (3.7) 



with 



TT "dcos9'JE'^ d COS 9'^dE'^ 



5,^iE,) = 5.ert(/^) + Cft(e,/i) + Cd(e,^,). (3.8) 

Here, d'^a^/ d cos 9'^dE'^ is the result of the numerical integration of the first term in (|3.4|) 
and the Born cross-section is: 



d COS 9',dE'^ S QU^, 



with the Born radiators 



S, = 2QI (3.10) 
2 



5: 



l-m)S'-M'Qf , (3.11) 



= ^-^Qf- (3.12) 
The hard part of the infrared divergent corrections is defined by Equation (1.4.8): 

1 r-^-T /■! /■27r E f O"^ \ 

5^^^^{e,E^) = - dE^E^ dcos9^ J^^^ = 2 In ^ In - 1 , 



TT Je J-1 JO e \ 



(3.13) 



where JF^^ is the Low factor (see e.g. Equation (1.4.5)): 



rlR _ 



The integrals used are given in the Appendix. 

The infrared divergent part of the soft photon contribution is: 

(e,/i) = - r dE^E^ dcos9^ d^^J^^^. (3.15) 

TT JO J-1 JO 
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With its calculation, one may follow Reference I until reaching Equation (1.4.26): 



Cft(e,/^) 



with 



V'^ + In ^] [{Q^ + 2m^)L^ - l] + S, + ^In^ + ^In^: 

/X '- /Pl i — Pl ^P2 J- — P2 



1 , v^ + Q: 
In 



2 

2 ' 
/ 



A. 



Qf(Q/ + W). 



(3.16) 



(3.17) 



(3.18) 



Calculating ( p.l6|) in the ultrarelativistic limit, with Pi and P2 being the velocities of the 
initial and final state electrons in the HERA frame, the infrared divergent part gets: 



cIR 

"^soft 



;6,/i) = 2(p'^ + \n-] (in^-lUln ^ 



+ 0$. 



(3.19) 



The function 5$ also reflects the dependence of the soft photon treatment on the reference 
system in which the corrections 61^^^ and S^^^ are calculated. We calculate 5$ numerically 
using the integral representation: 



1 da , I — da 
In- 



Pakl l + Pa 



with 



ka 



kia + ^2(1 — a), 



(3.20) 



(3.21) 



Pa 



Qfa{l - a) 



p^a + (pg) {I - a) + 2pePe cos 9^a{l - a) 



(3.22) 
(3.23) 
(3.24) 



E,a + E(l-a). 



(3.25) 



The infrared singularity is compensated by the contribution from the QED vertex 
correction (see e.g. Equation (1.4.34)): 



5vert(/i) = -2 ( + In - 1 ( In ^ - 1 1 - ^ In" ^ + ^ In ^ - 2 + — . (3.26) 



m 



Q 



2 Qi 



Qt 



TT 



In sum: 



2E^ 



1 



2\n--^\\r,:il,-l\-:-\n'^ + l\n%+\r.^^-2 + - + S^. 

2 m"^ 2 m'^ 



m 



TT 

"6 



(3.27) 
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In practice, we took = 10~^ GeV and verified stability against variations. With 
the Fortran program HECTOR, we further confirmed numerically with high precision that 
cross-section (|3.7|) agrees with that calculated in the semi-analytical approach provided 
no cut is applied. 



4 Higher-Order Initial-State Radiation 
in Leading-Logarithmic Approximation 

So far we considered the lowest-order cross-section. A calculation of higher-order photonic 
contributions is a non-trivial task since the cross-section is not completely inclusive. We 
will apply the leading-logarithmic approximation (LLA) which worked amazingly well for 
inclusive cross-sections (see, e.g., [|l]). 

We begin with the generic LLA formula for photonic initial state corrections. 



dz p{z) 



QED 



d COS d'^dE'^ Jo [d cos 6'^dE'^ 

The totally inclusive 0{a) density is: 



^ ^QEP ^ I. (4.1) 

^ d COS O'dEi I 



p(^) = ^i±i!('ln^-l). (4.2) 



2tt 1 — z \ m 
Only the initial electron energy scales, 

E,^E, = zE,. (4.3) 

The kinematical limits ( p.l6| ) and ( p.20| ) of the final state energies depend on the actual 
values of the angles and on Eg, Ep. Thus, we have to introduce a step function, 

e, = e{w^ - M^)e{Ml- M^), (4.4) 

in order to ensure the non-trivial kinematical conditions ([A.l|) and ( |A.2|) . Further, we note 
that the Jacobean is unity in (|4.1| ) since this differential cross-section is being calculated 
in terms of variables which do not scale with Ee — >■ zE^. 

How is the cross-section ( [4.1|) related to an experiment with photon tagging? Formally, 
one has to 'differentiate' ( ^.1| ) with respect to the visible photon energy. But we have to 
care about the definition of E'^^^. And this depends largely on the specific experimental 
set-up. There are two different cases to be considered: 

• Case (i) 

The tagging device is not in forward direction of the electron beam; 

• Case (ii) 

It is in forward direction of the electron beam. 
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The relevant object to be calculated is the three- fold differential cross-section 



dcos&JE'JE-^'^' ^ ■ ' 

where the visible photon energy is the energy deposited in the tagging device. 

We see two potential sources of inaccuracy which could make our calculations less 
reliable than usual applications of LLA often are: 

• For tagging in non-forward direction, the assumption of inclusiveness of the process 
in the LLA photon angle is fulfilled approximately while tagging in forward direction 
has to be treated with care (see below). 

• The photon, which is treated with the LLA method and that of the Born cross- 
section are identical particles but they are treated differently and their interferences 
are neglected. The treatment assumes them to be distinguishable. The resulting 
inaccuracy, though, is of next-to-leading order. 



4.1 Case (i): Non-Forward Tagging Device 

In Case (i), the visible photon energy is simply 



(4.6) 



The LLA describes the emission of collinear photons, here from the initial electron. These 
forward photons cannot hit the non-forward tagging device and are thus completely inclu- 
sive. Due to this, a direct application of ( |4.1j ) is justified as a reasonable approximation 
and the neglect of identity of the photons is expected to be not very influential. The 
cross-section with account of leptonic initial-state LLA corrections will be: 



i) 



d^a^j.. 



d cos ejE'dEZ 



d cos ejE'JEl 



+ 



d' 



a 



d COS e'JE'dEr' 



(4.7) 



with 



dcosOJE'JE:^'^ 



I d cos 6^d(p^ / dzp{z) 
JAn^ Jo 



d cos 6'^dE'^d cos O^dcp^dE. 



d cos 6'^dE'^d cos O^dip^dE^ 



^zEe 



4.2 Case (ii): Forward Tagging Device 

Case (ii) is realized in the HERA detectors. Here, the situation is much more involved. 
One could make the assumption that nearly every one of the initial state (LLA) photons 
is emitted collinearly, carries energy (1 — z)Ee and will hit the forward tagging device 
if its geometrical size is big enough. A simple estimate based on an analysis of the 
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denominator of the initial electron propagator after emission of the photon shows that 
the tagging device must have a geometric size fulfilling the condition 



AO, » 'f (4.9) 



-'e 



in order to collect most of these photons. For an opening angle A^^ = 0.5 mrad and 
Ee = 25 GeV this condition is fulfilled but not to such an extent that numerical estimates 
of the LLA corrections would be sufficiently well. In fact, the tagging device induces a cut 
on the photon angle and leads to an additional logarithmic dependence of the cross-section 
on it0. 

We have to distinguish several contributions. If both photons hit the tagging device, 

i.e.: 

< < A^^, (4.10) 
the visible photon energy has to be defined as the sum 

^vis,(iia) = (1 _ z)E^ + E^. (4.11) 

Next, if the photon of the lowest-order scattering process does not hit the photon tagger, 
i.e.: 

> Ae^, (4.12) 

only the coUinear LLA photon contributes to 

= (1 - ;z)Ee. (4.13) 

Finally, it may happen that the LLA photon does not hit the tagger, but the lowest-order 
photon does: 

0<ej< Ae^, (4.14) 

and the visible photon energy is 



All three cases contribute to the cross-section and have to be summed up with the 
lowest-order cross-section: 



dcosO'JE'JE^'^ dcosO'JE'JE^'^ dcosO'JE'JE:;^'' 

d^(^iiib) d^a^-i^^ 
dcosO'JE'JE:^^' dcose'JE'JE:;^'' ^ ' ' 

^ LLA corrections with photon tagging have been discussed in where a bremsstrahlung cross- 
section, reduced by tagged events, was estimated. The aim was to reduce the radiative corrections from 
hard photon emission. Our modifications of the density function by photon cuts are quite similar to the 
treatment in [^. 
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We now have to extract E^^^ from ( f4.1| ). This will be done in the next three subsections 



with use of identity 5-function 


integrations: 




1 = 


/ dEf5 


"Sf - 


-(1 


1 = 


f dEfd 




-(1 


1 = 


f dEfd 


"^vis _ 


- E^ 


The three cases will be treated 


separately. 





-(l-z)E, 



(4.17) 
(4.18) 
(4.19) 



4.2.1 Case (iia): Both photons hit the forward tagging device 

When two photons hit the forward tagging device, the cross-section gets: 



dcosO'JE'JE:^'^ 



dz pa{z) 



d^ 



a 



X 



QED 



d COS 6'^dE'^ 



dcosO-^ 



d'^cr^^^ 

Q _ QED 



dip^ / dzpa{z) 



d cos O'^dE'^ 



;i - z)Ee - E, 











X 



+ 



dE^ 



d cos 9'^dE'^d cos O^dip^dE^ 



-0(£/'y — E^ 



d^a 



An d cos 9 ^dE^ 



dE^ 







Ee^zEe 

5, 



(i (Tbr 



d COS 6'^dE'^d cos 6^d(p^dE^ 



-Q(^Ey — E^ 



An a cos o^aE^ 



>6 



(1 - z)E, - E^ 



As mentioned, the density Pa{z) will deviate from the inclusive density 
angular cut: 



(4.20) 
due to the 



Pa{z) 



271 1 — z \ m? 



In- 



- 1 



(4.21) 



After performing the integral over z and some trivial re-orderings, we arrive at 



dcos&JE'JE:^'^ 



A 9, 



7 







1 r2-iT /■mm{£;:!^'=,£;iji="'} 

2 dcosO^i I dcoscp^ I , ^ dE, 



Pa{z) 

E. 



2 = 1- 



^7 --^7 



d Tbrem 



7 



dcos O^dE^dcos 6^d(p^dE^ 



Ee^zEe 



d^cr^r 



dcos O^dE^dcos O^dip^dE^ 
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+ 









47r E, 


civis 

z ■ -1-^ 


d cos 9'^dE'^ 



d COS 6' dE' 



0(^min — Zq) 



(4.22) 



Everything is straightforward in the above derivation. The only dehcate point is the 
determination of the boundaries for the integration over E^, and, related to this, the 
allowed range for E^^^. This will be explained in some detail now with use of Figure |^. 
The integration over E^ proceeds along one of the dotted straight lines with defining 
equation (|4.17| ): 



E^ 



Ee]+zE, 



(4.23) 



The absolute maximum of EZ^^ is reached for z = 1, and the absolute minimum for: 



E^ 
E,, ■ 



(4.24) 



where E^ vanishes. For the hard photon integral in (|4.22|) this minimum has to be replaced 
by E^ as discussed above. The discussion of limits would be finished here if there wouldn't 



be the scaled kinematical upper limit for E^, 
limit depends on several variables: z, E^, EL 



TTimax 

Z-, E^, Ep, Eg, ^g. 



which also has to be respected. 



This 
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From ( |2.2(J| ) one may see that the minimum value of -E™^^, which is zero, is independent 



of the actual values of the photonic angular variables. It depends, besides on Ep, on zEe 
and external variables E'^,6'^ and is reached for fixed E'^,9'^ at some value zq which may 
be easily derived from ( |A.4| ) (with neglected terms of order 0{m^)): 



zq = zo{Ee, Ep, Eg, 



2ki{pi - k2) 



(4.25) 



may be bigger or smaller than 2;min in dependence on 
mm, it is of no influence on the integration region and the 



As is visualized in Figure ^, 
E^^^. If zq is smaller than 2;, 
lower integration boundary is E^. In the other case, one has to study in detail where the 
necessary condition E'^^^ > i?™ is fulfilled. Some value zi, corresponding to some value 
^min jjiay be found where this begins to be fulfilled and the lower integration boundary 
becomes E. 



"■mm. 
"7 • 



^7 



E^''^{ziEe, Ep] Eg, (p^, 6'^] 



1 - zi)E, 



(4.26) 



Therefore, d cr^^ will not contribute to ( |4.22| ) 



We further have to study the fan of curves defining the value of E™*^^, starting from 
Zq and depending on the parameters 6^, ip^ as a function of z at fixed E'^,9'^. Figure || 
shows that, for HERA kinematics, the curves are nearly straight lines. If their curvature 
is strong enough, they may cross the straight dotted integration path from above, in 
principle independent of the relation between zq and ^min- If this happens the upper 
integration limit in ( [4.22| ) is not E™ but E™*™. The figure suggests that for the case of a 
forward tagging device (with cos^^ ~ 1) the relevant solid curve is nearly parallel to the 
dotted line. Then, for reasonable choices of E^'^ no problems should occur from details 
of the kinematics. 

To conclude the discussion, one has to integrate the scaled cross-section over all values 
of E^ on one of the dotted lines which are located below the envelop of the fan of curves 



E^ 



The unsealed cross-sections are integrated from E^ to -E™. In applications, we 



restrict ourselves to kinematics with zq < z, 

stops if Zj^ia < Zq. 



mm and 



> 



The FORTRAN program 



4.2.2 Case (iib): The LLA photon hits the forward tagging device 

If only the collinear initial-state LLA photon hits the forward tagging device, the cross- 
section becomes: 



dcose'JE'JE:^^^ 




QED 



d cos d'AE' 



'1 - z)E, 



2tt 



dip^ I dz pb{z) 
Jo 

d Cbrem 



d COS 9'^dE'^d cos 6^d(p^dE. 



■Q{E^ — E^ 
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+ 



d'^a. 



d^abr 



d cos 6'^dE'^d cos O^dip^dE. 



-Q(^E^ — E^ 



>6 



- (1 - z)E, 



An d cos O'^dE'^ 

The density function for the LLA correction is the same as that for Case (iia): 

Pb{z) = pa{z). 



(4.27) 



(4.28) 



The subsequent formal steps differ from Case (iia) since E^ is not involved in the 6- 
function. The limits of variation of E^ are affected only by scaling: 



d'a] 



(iib) 



dcosO'dE'dEl 



Pa{z) 

E. 



(icos6'^ / dcosip-y 

1 JO 



dE^ 



d'> 



dcos O^dE^dcos O^dip^dE^ 

d 0"brem 



dcos O^dE^dcos 6^d(p^dE^ 



+ 



An 



d'^a. 



d cos 6'dEl 



-Ee *^min-^e 



d cos 6' dEl 



(4.29) 



The integration boundaries may be inspected in Figure |^. One should notice that the 
upper integration limit is again different for the scaled and the unsealed cross-section 
integration. In (|4.29|) , the unsealed integration limits are used, however. In the figure, 
they define the solid straight line ranging from the point (z, E^) = (1, E^) upwards until 
it meets the actual value of E^^^; the latter depends on the kinematics. The scaled 
limits are dependent on z and define a dotted straight line at fixed Zmin- These limits are 
automatically respected by means of the factor Qz^i^ ( [4.29|) . 

Further, we introduced a new variable, E^^^, in the definition of the upper integration 
bound. If somebody wants to impose a general cut on the maximum allowed photon 
energy, it has to be done here and only here. 



4.2.3 Case (iic): The LLA photon does not hit the forward tagging device 



This case is very similar to Case (i). The correction may be read off from (|4.8| ) with a 
modified density pc due to the cut on the angular integration for the LLA photon and 
with performing the angular integration for the Born photon as in Case (iia): 



{iic) 



dcosOJE'JE:^'^ 



dcosOy 



2n 



dif^ / dzpc{z) 



d^(Thr 



d cos 9'^dE'^d cos O^dip^dE^ 



*zEe 
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d C"brcm 1 

d cos cos O^dip^dE^ J ' 



(4.30) 



with 



5 Numerical Results and Discussion 

Numerical results are shown for the HI forward tagging device, which is approximated 
here by a circle located symmetrically around the beam axis with A6'^ = 0.5 mrad. 
The lowest-order cross-section has been shown in Section ^. The estimated radiative 
corrections (5, 

d^a'^^'/dxdydEZ'' 

' ^ ^ X 100%, (5.32) 



d^a\,T-em/dxdydE, 



VIS 



are shown in Figures ^ and |^. The d^a^^^'^ is defined in ( p.34| ) and d'^cr™' in ( |4.7| ) or 
(HI)- 




^min ^min 



Figure 5: The physical region (z, Ej) for Case (iib): only the LLA photon hit the forward 
tagging device 

For the figures, we consider Case (ii) (forward tagging) and assume i?"*^ = 5 GeV and 
E'l^^ = 10 GeV. The maximal allowed photon energy is assumed not to be restricted; 
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this corresponds to no photon observation outside the tagging device. As was mentioned 
already in the discussion of the Born cross-section, there is a minimal value Uc of y which 
depends on the kinematics Otherwise, the resulting corrections are similar 

to those known from the inclusive measurements. The Born cross-section dies out for 
X > 10~^, thus limiting the values at HERA to be not bigger than several GeV^. For 
nearly all this kinematical range, the corrections stay well in between about ±20%, mostly 
even within and 10%. Larger, negative QED corrections arise in the soft photon corner 
(at small y) and large positive ones for hard radiation at high y. The biggest, negative 
corrections are obtained at larger values of x where the Born rates go down rapidly. The 
soft photon corrections are weakened by higher-order contributions. These may be taken 
into account by performing soft photon exponentiation. An approximation could be to do 
this for either the lowest-order bremsstrahlung cross-section or the LLA density function 
separately. An explicit higher-order calculation is needed for a sufficiently correct answer. 
The numbers in the figures were produced without soft photon exponentiation. 

In this article, we studied deep inelastic cross-sections at HERA energies with photon 
tagging. The lowest-order radiative process is calculated in ultra-relativistic approxima- 
tion with no additional restricting assumptions. The next-order corrections have been 
estimated with a leading-logarithmic approximation neglecting the interferences between 
diagrams with different topology concerning the two photons. The resulting error is ex- 
pected to be of next-to-leading order and could be non-negligible for a forward tagging 
device. Our calculations are performed strictly in order 0{a). For soft photon emission, 
we know to overestimate the large, negative corrections. An exponentiation of soft photon 
corrections is recommended here if one wants to improve the numerical results. Again, 
problems due to interference effects remain unsolved so far. For a forward tagger, the 
problem is not completely inclusive. Logarithmic corrections due to the finite size of the 
tagger are taken into account. 

Since the calculated corrections are not too big compared to the typical size of photonic 
corrections in deep inelastic scattering one may expect that they represent a quite good 
approximation of a complete calculation of the order 0{a) photonic corrections. 
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Figure 6: Photonic correction 5 = d^a'''' / d^a^,^^ x 100% for E:^^ = 5 GeV and M. 
mrad for a forward tagging device 
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A Kinematics and Phase Space 



The angular variables 9'^, 9^, Lp^ are unrestricted. The lower limits of E'^ and E^ are m 
and zero, respectively. 

The upper limits of the two energies depend on the order of integration. We calculate 
the cross-sections in terms of the electron variables E'^,9'^. Therefore, we call them external 
variables and may derive the corresponding boundary from the condition 

= (E, + Ep-Ely- (pf - {p, -ppy + 2p2pe - Pp)cos 9'^ 

> M\ (A.l) 



Solving this inequality for E'^ yields ( |2.16 



The two-dimensional region of variation of cos 9^ and E^ is shown in Figure ^. Photonic 
variables are considered to be varying at fixed values of electron variables and will be 
called internal variables. The limit for E^ follows from a condition which relates /c2 and 
k properly. It is: 

Ml = {pi + ki-k2- kf > M\ (A.2) 
In terms of natural variables: 

Ml = (A.3) 



-2E^< 



PgCos 9'^ — (pe — pp) COS 9^ + PgSin ^gSin 6'^cos — E'^ + E^ + Ep> , 



from which we get 

2 [Ee + Ep- E'^ + (pg cos 9'^- \pe- Pp]) cos 9^ + p'^ sin 9'^ sin 9^ cos cp^ 



which in fact is ( p.20| ). 

In the ultra-relativistic approximation in m, (|2.2CI|) reduces for 6*^ = to 

™ ^E^Ep + (M^ - M')/2 - E^E'^jl - cos^l) - EpE'^jl + cos^l) 

^7 = ^ — — z:z^\ ' (A.5) 



2Ep-E'^{l- COS 9'^'j 



and for 6*^ = vr to 



2E,Ep + (M^ - M')/2 - EXa ' cos^l) - EpE'^jl + cos^l) 
^ 2E,-E;(1 + cos^;) ' ^ ' 

These limiting cases are of use for an understanding of the plots of kinematical boundaries 
of variation of the internal variables cos^^, (p^, E^ which are shown in Figure 

Next we have to discuss the phase space parameterization. We aim at a check on 
the correctness of the kinematical boundaries since the angular bounds have not been 
explicitly derived. 
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COSl^j 



Figure 8: Kinematical boundaries for cos 9^ and 

In Eq. (1.3.9) we made use of a parameterization of the phase space in terms of the 
invariants yi,Qf,yh, Ql, together with one additional variable Z2 (or zi): 



dT 



—-dyi dQi dyh dQf, 



(A.7) 



From explicit expressions for yi and Qf in terms of natural variables it is easy to derive 
the Jacobean 



dyidQf = 2p^dE'^dcos6'^. 



(A.8) 



From (|A.7| ) and (|A.8| ) we get the doubly-differential phase space in terms of cos9'^,E'^ 
from a three-fold integration over invariants with the aid of Eqs. (I.D.2), (I.C.2), (I.C.3): 



dT 



p^dE^dcosU^ 2 Jm^ •'Qlmin •'""^^^ t^vR 



(A.9) 



where we also used (1.3.7): 



dyhdQl = -dMldQl. 



(A.IO) 



On the other hand, the phase space ( |A.9|) in terms of natural variables was given in ( p.23|) : 



dV 



IT 



p'^dE'^dcosO'^ 2 Jo 



2lT 



dif^ 



dcos 6-, 



^max 



(A.ll) 
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a) cosi5;= 1 , E;= 1 5 b) cosi^;=0, E;=25 




c) cosi?;=-i, e;=400 



Figure 9: Kinematical boundaries for cos9^,ip^, at given values of E'^ and cos 9'^ 
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We checked numerically that the phase space ( [A .91 ) agrees within computer precision (14 
digits) with the new expression (|2.23|) if one uses in (|A.9|) the W'^ introduced in (|A.1|) . 
Finally, we have to collect the integrals for the calculation of S^^^: 



[A] 



TT 



d cos 9^ 



2-K 



d(p-yA. 



(A.12) 



For sufficiently small E^ there is no dependence of the angular integration limits on E^. 
In this case, the integrations are straightforward and the results are very simple: 



2 ' 

m 




^2 





m 



^2 
^2 



Z1Z2 



In^, 
e 

2L,,ln^, 
e 



(A. 13) 
(A. 14) 



and Ljn is defined in ( |3.171) . The above two integrals are the equivalent of the tables of 
integrals in appendix (I.D.3) for the present calculations. 
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